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PREFACE 


This report provides coastal engineers algorithms for a number of hand- 
held calculator programs for coastal engineering, primarily in the area of 
wave transformations and wave generation. These algorithms were developed 
under the coastal engineering research program of the U.S. Army Coastal 
Engineering Research Center (CERC). 
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CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (SI) UNITS OF MEASUREMENT 


U.S. customary units of measurement used in this report can be converted to 
metric (SI) units as follows: 


Multiply by To obtain 
ims Ln 25.4 millimeters iy es 
2-54 centimeters 
square inches 6.452 square centimeters 
cubic inches 16.39 cubic centimeters 
feet 30.48 centimeters 
0.3048 meters 
square feet 0.0929 square meters 
cubic feet 0.0283 cubic meters 
yards 0.9144 meters 
square yards 0.836 Square meters 
cubic yards 0.7646 cubic meters 
miles 1.6093 kilometers 
square miles 259.0 hectares 
knots 1.852 kilometers per hour 
acres 0.4047 hectares 
foot-pounds 1.3558 newton meters 
miblele bates 1.0197 x 1072 kilograms per square centimeter 
ounces 28.35 grams 
pounds 453.6 grams 
0.4536 kilograms 
ton, long 1.0160 metric tons 
ton, short 0.9072 metric tons 
degrees (angle) 0.01745 radians 
Fahrenheit degrees 5/9 Celsius degrees or Kelvins! 


1To obtain Celsius (C) temperature readings from Fahrenheit (F) readings, 
use formula: C = (5/9) (F -32). 
To obtain Kelvin (K) readings, use formula: K = (5/9) (F -32) + 273.15. 
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HAND-HELD CALCULATOR ALGORITHMS FOR COASTAL ENGINEERING 


by 
Todd L. Walton, Jvr., Willtam A. Birkemeter, and J. Richard Weggel 


I. INTRODUCTION 


The advent of the hand-held programable calculator has led to the develop- 
ment of numerous programs in various fields of engineering and science. 
Coastal engineering is no exception. This report provides algorithms for a 
number of calculator programs useful in performing coastal engineering calcu- 
lations, primarily in the area of wave transformations and wave generation. 


There are basically two types of hand-held programable calculators: those 
that use algebraic logic, such as Texas Instruments, Algebraic Operating System 
(AOS) notation, and those that use Reverse Polish Notation (RPN), such as 
Hewlett-Packard. For the six programs presented herein, versions in each type 
of notation are provided--AOS designated by "A" and RPN designated by "R." 

Each program is documented, the assumptions are briefly described, and refer- 
ences to more detailed presentations of the theory are given. 


Since subtle differences exist among various brands of calculators and 
among various models by the same manufacturer, the programs presented may need 
to be modified to adapt them to other brands or models. These programs were 
prepared and checked, using the Texas Instruments TI-59 (AOS) and the Hewlett- 
Packard HP67 (RPN). Users should refer to the operating manuals supplied with 
individual calculators. Algorithms may be set up specifically to meet some 
limiting conditions. For example, some programable calculators have a limiting 
restriction on the number of program steps accepted so a program must be re- 
vised to meet the limitation. The Texas Instruments TI-57, for example, 
accepts only 50 program steps. Algorithms may be adapted to meet this limita- 
tion by requiring prior calculation of input parameters either manually or by 
another program. 


Each of the AOS programs incorporates TI-59 compatible print routines which 
print and label all input and output parameters. The user only has to enter 
the input parameters and the results are automatically computed and printed. 
Since the printing routines increase program length by as much as 25 percent, 
use of a magnetic card for permanent program storage is recommended. A11 
print steps are marked with asterisks and need not be entered if printing is 
not desired. 


II. PROGRAMS 


Six programs (100, 101, 102, 103, 104, and 105) are presented in this 
report. Program 100, a simple program that computes linear wave theory wave— 
length for a given depth, is designed to be used as the basis for any program 
that requires wavelength; in fact, it has been incorporated into programs 101, 
102, and 105. 


Program 101 is another basic. program which computes not only wavelength 
but also a number of other linear wave theory parameters. This program forms 
the basis for program 102 and can be adapted to other programs as well. 


Program 102 computes linear wave parameters and breaking wave height and 
direction based on nearshore or deepwater wave information. Program 103 can 
be used to forecast wave height and period in shallow water. Program 104 and 
105 address wave conditions at structures—-program 104 predicts the depth- 
limited design breaking wave height at a structure; 105 uses Fuch's equation 
to predict wave transmission over a thin barrier. 


Each program allows either English or metric input and output (the AOS 
programs print an "M'' when metric is selected). Program listings are anno- 
tated, making it possible to follow the logic of the algorithm and to make 
modifications if desired. 


There are undoubtedly many calculator programs not included here that have 
been developed on coastal engineering subjects. Practicing engineers who would 
like to disseminate such programs (in either AOS or RPN) to other users are 
encouraged to submit them to the Coastal Engineering Research Center (CERC). 

If the response is great enough, additional reports presenting the programs 
will be prepared. Comments, programs, or suggestions for programs should be 
sent to: 


Commander and Director 

US Army Coastal Engineering Research Center 
ATTN: Evaluation Branch 

Kingman Building 

Fort Belvoir, VA 22060 


These programs and future programs will generally correspond to the following 
numbering scheme: 


Miscellaneous 0-99 
Waves and currents 100-299 
Inlets 300-499 
Beaches 500-699 
Geology 700-899 
Structures 900-1099 


In general, the documentation of programs submitted should be in a format 
paralleling that of the programs presented in this report. A blank set of 
forms which can be reproduced is included in the Appendix. 


Program Description 


Program Title 100R Linear Wave Theory Wavelength (RPN logic) 
Name J.R. Weggel Date 12/80 
Coastal Engineering Research Center 

Kingman Building “ee Zip Code 
Fort Belvoir e 


Address 
City 22060 
Program Description, Equations, Variables, etc. 
This algorithm takes deepwater wavelength as input and using the depth 
at a given site iterates to obtain wavelength by linear wave theory. 
Algorithm uses English or metric system of units. 


REFERENCE 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, 
Shore Proteetton Manual, 3d ed., Vol. I, Eq. (2-4), Stock No. 008- 
002-00113-1, U.S. Government Printing Office, Washington, D.C., Oe 


Operating Limits and Warnings 
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User Instructions 


100R Linear Wave Theory Wavelength (RPN logic) 
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Program Description 


Program Tite 100A Linear Wave Theory Wavelength (AOS logic) 

Name J.R. Weggel Date 12/80 
adress Coastal Engineering Research Center 

Kingman Building 


t : 
Gu Fort Belvoir, 


State Zip Code 


Virginia 22060 


Program Description, Equations, Variables, etc. 


This algorithm takes deepwater wavelength as input and using the depth 
at a given site iterates to obtain wavelength by linear wave theory. 
Algorithm uses English or metric system of units. 


REFERENCE 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 
Protection Manual, 3d ed., Vol. I, Eq. (2-4), Stock No. 008-022-00113-1, 
U.S. Government Printing Office, Washington, D.C., 1977. 


Operating Limits and Warnings 


100A-1 


User Instructions 


Program Title: 100A Linear Wave Theory Wavelength (AOS logic) 
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Program Description 


Program Title 101R Calculation of Wave Parameters from Linear Theory(RPN opie 
Name T.L. Walton, Jr. Date 12/80 
Coastal Engineering Research Center 
Kingman Building 

: Zip Code 
Fort Belvoir irgind i 22060 


Address 
City 


Program Description, Equations, Variables, etc. 


This program calculates the product of the wave number and depth, kd, the 
ratio of group wave speed to wave celerity, n = 0.5 (1+2kd/sinh 2kd), the 
group wave speed, Cg, the shoaling coefficient, Ks, the refraction coeffi- 
cient, Ky, horizontal orbital velocity, u, and vertical orbital velocity, w. 


Program input includes wave period, T, deepwater wave angle, 4), deepwater 
wave height, Ho, wave phase angle, 8, depth of water, d, at which results 
are desired, and depth from surface, z, at which velocities are calculated. 
This program assumes straight and parallel offshore bottom contours for 
assumption of Snell's law of refraction. Algorithm uses English or metric 
system of units. 


REFERENCES 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 
Proteetton Manual, 3d ed., Vol. I, Ch. 2, Stock No. 008-022-00113-1, 
U.S. Government Printing Office, Washington, D.C., 1977. 


Operating Limits and Warnings This program can also be used for wave input H, T, 
and a, from transitional water depths (as in the case when wave gage infor- 
mation is provided but the gage is not in deep water). In this instance, 
some program steps must be deleted as noted in the program listing, and the 
Kr computed in the program is not the refraction coefficient for depth, d, 
but a modifted refraction coefficient from the gage site to water depth, d. 
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User Instructions 


101R Calculation of Wave Parameters from Linear Theory 
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User Instructions 


_1Q1R Calculation of Wave Parameters from Linear Theory 
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Program Description 


Program Title 101A Calculation of Wave Parameters from Linear Theory(AOS logic) 
Name W.A. Birkemeier Date 5/81 
Address CERC Field Research Facility 


Kitty Hawk, State North Carolina Zip Code 27949 


City 
Program Description, Equations, Variables, etc. 

This program calculates the product of the wave number and depth, kd, the 
ratio of group wave speed to wave celerity, n = 0.5 (1+2kd/sinh 2kd), the 
group wave speed, Cg, the shoaling coefficient, Kg, the refraction coeffi- 
cient, Ky, horizontal orbital velocity, u, and vertical orbital velocity, w. 


Program input includes wave period, T, deepwater wave angle, 09, deepwater 
wave height, Ho, wave phase angle, 9, depth of water, d, at which results 
are desired, and depth from surface, z, at which velocities are calculated. 
This program assumes straight and parallel offshore bottom contours for 
assumption of Snell's law of refraction. Algorithm uses English or metric 
system of units. 


REFERENCE 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 
Protection Manual, 3d ed., Vol. I, Ch. 2, Stock No. 008-022-00113-1, 
U.S. Government Printing Office, Washington, D.C., 1977. 


Operating Limits and Warnings This program can also be used for wave input H, T, and 
a, from transitional water depths (as in the case when wave gage informa- 
tion is provided but the gage is not in deep water). In this instance, some 
program steps must be deleted as noted in the program listing, and the K; 
computed in the program is not the refraction coefficient for depth, d, but 
a modified refraction coefficient from the gage site to water depth, d. 
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User Instructions 


Program Title: 1Q1A Calculation of Wave Parameters from Linear Theory 
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AREER S SCaA 
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User Instructions 


Program Title: 101A Calculation of Wave Parameters from Linear Theory 


INPUT OUTPUT 
ITH PRinTeR. WyU, AND WARE PRINTED _IN Sieg) 
APPROPRIATE UN ITS (VELOCITIES ARE Pi ecome een 


° 
> 
2 
< 
z 
a 


“ 
0 
m |B Ie 
r. 
i 


{3 CTKERWISE? 


READ H IN DISPLAY (CIF USING YU A essen 


H 


READ U-IN DISPLAY AND PeRocEED TO step F) 
PRESS R/S AND READ U IN DISPLAY 


u 


PRESS: R/S AND READ WIN DisPLA WwW 


To REPEAT, START AT STEP Foe NEw © oR 


TH TeP A cR B® AND TH 


JUCGOUUL UCU 0000s 000000 
JUUUCUUUUUUU UCU 020 000000 
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User Instructions 


Program Title: 101A Calculation of Wave Parameters from Linear Theory 


T= &8sec 
d= Soft 


Ho = IBft 


NTER Sy PRESS © 


PRESS E To ComPuTe 


CALCULATED VALUES ARE PRINTED GR 


| | PRESS R/S AND ReaD Cq = 24,625 ft/sec 
| | Press R/S pup READ Ka= 0.4I2 
| | press R/S AnD READ y= 0.095 


EnTER Zz 
ENTER C© peess C' 


CALCULATED VALUES ARE PR)NTED OR: 


oma 
= @ e 


hou 


hw htt 


QOPI ON 
CT Boy bp 


MMe wo 


. 
AM hae 
DOW hw om 
FP hOM OO 
Ou Aa 


AUNOU 
DOGO Saas 


HPUMHee 


apa 
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Program Listing 


STEP KEY ENTRY KEY CODE COMMENTS STEP KEY ENTRY KEY CODE COMMENTS 
nog Fe LE ey SI Nae SIN %_ IN DISPLAY 
(ami = shal fe as gta (Wale 
oz 42 $To TR, O56 38 SIN Magara ane? 
ate) re oS? 42 STO Oe 
fae Meas * nes OF! wr sinihic™) 
ioe Mr oz % "+" pet CODE (aye) eee TRY 
no6 63 OF * 060 23 LH 
oo7 G4 O4 »* May  7R s 
noe 43 RCL «* M62 43 RCL 
ooo of oe «x Gee Wr Wr 
nig «69 OF x 64 94 4¢7- 

Hie Os Gece es 65 22 TH 
pees a Rasy n6e 23 LHe 
13 76 LEL iar a ay 
aheh shee fs} OS Be + 
O1S 42 STO dR Hes) We 2 
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nig oo6 4 x wre Fe (Lei 
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q20 m2 aa O74 42 STO Conte 
Da AS Re Gres OF) Or cosh (x) 
Me fi ml co of6e 22 IH 
Mee teas) (al? € rm Ze line 
(22 @3 “Ro. Yabuc ics ese 
G25 65 x OF3 43 RCL 
026 69 Oo) ly wry 
ay (Say Re G81 94 +7- 
Os He 2 os2 22 IH 
29 a5 = uss 23 LH 
20 42 570 gad > Ru 034 33 = 
iol el aes. 55 = 
ie Ole CE Ue te a ge 02 2 
33 76 LBL Wey Ss 0S 
geet. ‘lay os8 92 RTH 
Osseo oH 7 Ro o69 76 LBL 
036 oO OO o90 i4 OD 
ey i % g a 4 
Of OF 3 «* “ANGO” PRINT Code O93 O8 & ae poate 
O40 01 1 % 034 o1 1 
Od a o95 42 STO 
O42 02 2 « Os 18 ie DELETE IF 
43 of 0 «* Wer GE ~E ose UNWANTED 
O44 O01 1 x USS OC %K 
045 69 OF « age ee ipo Gini?! 
O46 04 O4 x 100 04 o4 #* Rix 
O47 43 RCL «x LO wes) CF. * 
042 00 OG + Oe | Oe (Oa eee ee 
049 69 OP % = PRINT Xo 103 61 GTO INDICATE METEIC 
O50 06 06 * 104 4 (apie 
O51 38 SIH 105 76 LBL 
as2 42 STO sin Y Ro NOS We 
153 \) oo t 
101A-5 


% - DELETE IF PRINTER 13 NOT AVAILABLE 
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STEP KEY ENTRY KEY CODE COMMENTS STEP KEY ENTRY KEY CODE COMMENTS 


3- 32.2 F b/sec® 


ITERATE UNTIL 
Luew-Locp st 
WHERE t IS VALUE 
IN tC REGISTER 


Oe 
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OMe OVO yO bh hae 
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M-DELETE 1F PRINTER IS NOT AVAILAGLE 
R/S- INSERT IF PRINTER IS NOT AVAILABLE 
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STEP KEV ENTRY KEV CODE COMMENTS STEP KEV ENTRY KEY CODE COMMENTS 


4 271 at st hal 
= BEE ar => 
219 42 370 n>Rio zr4 03 O3 
220) 10) #10 27s) | Sie) 
Bet ays) (oP PRINTN #76 43 ROL 
Z o6 Aya (oe) we) 
= 3a erie | Se 
224 1 + 279 33 KE 
225 5 « “eg. PRINT CODE 230 75 = 
22 baie es oe ZS Oil 
By OR BG 282 95 = 
22s 69 OF «* 283 394 ¢+-- 2 
299 04 O4 » 284 25 1-7 Coss 
230 43 RCL 25) Pos es 
ee inet (mel DES «Gel & 
Se SS) a 267 O1 1 
233 43 RCL ete Le eS 
244 02 2 239 43 REL 
Bei SG) 2 290 OO Og 
236 69 OP * prt Cg A 33 HE 
227 06 6 2 ea = 
WS See. et se ACE) Bes Ire 
239 65 x 294 34 THR 
240 02 2 % e295 42 570 
aol Me 8B ft ooo eS O eee wl yaa 
Som O68 2  G Be Mantz 297. 69 UPL aase Wo eu 
243 06 65 * 298 06 06 « 6 
e44 69 OP * 299 98 ADDY 
245 04 o4-* 300 91 Res ae 
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247 2 fe eye ey Tale REQUIRED IF CRRITAL 
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Para) ey abe ce We. 2 “ to 7 Re 
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252 04 4 e307 oo 4 « 
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259 69 OP * 214 69 OF *« PRINT Ho 
S30) (3 dia re Seow Xe 315 06 6 * 
aa We Be e 316 43 ECL 
202. 06 & 8 m0 O gily iti alia 
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260) 65 ORD x e 320 O11 O1 
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263 00 0 O eyes} (0) (me) 
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* -DELETE IF PRINTER ID NOT AVAILABLE 101A-7 
R/S— INSERT IF PRINTER ISNO6T AVAILABLE 


@ — DELETE IF Hy) NOT Ho IS KNOWN 


27 


STEP KEY ENTRY KEY CODE COMMENTS STEP KEY ENTRY KEY CODE COMMENTS 
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K — DELETE IF PRINTER 15 UNAVAILABLE 101A-8 
R/S — INSERT IF PRINTER 15 UNAVAILARLE 
© — DELETE IF H, NOT Ho 1S KNOWN 
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Program Description 


102R Linear Wave Approximation to Breaking Wave Height and 
Program Title Breaking Wave Angle (RPN logic) 


Name T.L. Walton, Jr. Date 10/80 


Address Coastal eee Research Center 
Kingman Building Zip Code 
ey Fort Belvoir, Stale Virginia Me 22060 


Program Description, Equations, Variables, etc. 
This program calculates breaking wave height, Hp, and breaking wave angle, 
db, using linear wave theory approximations combined with the shallow-water 
breaking assumption. Input parameters are wave height, H, wave period, T, 
wave angle, a, and the water depth, d, where the preceding three variables 
are measured. An additional input parameter is nearshore beach slope, m. 
The ratio of the breaking wave height to the water depth at breaking is pre- 
dicted using the equation 


m 
K = Hp)/dp = 1.16 TEL) 


from Singamsetti and Wind (1980), where dj} is the water depth at breaking, H! 
the deepwater wave height, and Lo the deepwater wavelength. This solution 
requires the assumption of straight and parallel offshore bottom contours for 
the application of Snell's law of refraction. Input wave parameters H, T, 
and a can be in any depth of water, d. Algorithm uses English or metric sys 
tem of units. The development of the equation is derived on the attached 
solution sheet. 


0.22 


REFERENCES 


SINGAMSETTI, S.R., and WIND, H.G., "Characteristics of Shoaling and Breaking 
Periodic Waves Normally Incident to Plane Beaches of Constant Slope," 
Report No. M1371, Toegepast Onberzoek Waterstaat, July 1980. 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 
Protectton Manual, 3d ed., Vol. I, Ch. 2, Stock No. 008-022-00113-1, 
U.S. Government Printing Office, Washington, D.C., 1977. 


Operating Limits and Warnings 


102R-1 
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SOLUTION SHEET FOR PROGRAM 102R 
Development of the equation: 
From conservations of energy 
2 H;2 
- Cy cos a = 1 Cei COS af @) 
where the subscript i indicates incident wave parameters. 


If left hand side of above equation represents conditions at breaking then 


Cg = € = Ch = Yedp = yshb/k (2) 
Hp 
where Cee (3) 
b 
Now assume K = 1.16 pe Nae (4) 
Bo Gus) 


where Hj is unrefracted deepwater wave height. 


Using (1), (2), (3), and (4) it can be found 
1/2 2/5 
Hp -{() / H,? Ceti cos a: / (5) 
From Snell's law of refraction 


sin op) sin ot 
Qi Gk 


therefore, 


1 
oom ache ° 


102R-2 
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User Instructions 


102R Linear Approximation to Breaking Wave Height 


and Breaking Wave Angle 


Program Title: 


OUTPUT 


*~DATAJUNITS 


INPUT 
OATAJUNITS 
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102R-3 
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iii ome Listing 


KEY ENTRY KEY — 


STEP 


COMMENTS 
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COMMENTS 
SOBBQOVAIWE 
Cosi Cc ) 


KEY CODE 
ieSio2 


i) 
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KEY ENTRY 
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Program Description 


ae 102A Linear Wave Approximation to Breaking Wave Height and 
rogram Title Breaking Wave Angle (AOS logic) 


Date 5/81 


PTO W.A. Birkemeier 
Address CERC Field Research Facility 
City Kitty Hawk, State North Carolina ZipCode 27949 


Program Description, Equations, Variables, ete. 
This program calculates breaking wave height, Hp, and pene wave angle, 
Ap, using linear wave theory approximations combined with the shallow-water 
breaking assumption. Input parameters are wave height, H, wave period, T, 
wave angle, a, and the water depth, d, where the preceding three variables 
are measured. An additional input parameter is nearshore beach slope, m. 
The ratio of the breaking wave height to the water depth at breaking is pre- 


dicted using the equation 
/ ane m 0.22 
K = Hp/dp = 1. a) 
b/ ap VHS/Lo 


from Singamsetti and Wind (1980), where dp is the water depth at breaking, 
H} the deepwater wave height, and Lo the deepwater wavelength. This solu- 
tion requires the assumption of straight and parallel offshore bottom con- 
tours for the application of Snell's law of refraction. Input wave para- 
meters H, T, and a can be in any depth of water, d. Algorithm uses English 
or metric system of units. The development of the equation is derived on 
the solution sheet included with program 102R. 


REFERENCES 


SINGAMSETTI, S.R., and WIND, H.G., "Characteristics of Shoaling and Breaking 
Periodic Waves Normally Incident to Plane Beaches on Constant Slope," 
Report No. M1371, Toegepast Onberzoek Waterstaat, July 1980. 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 
Protection Manual, 3d ed., Vol. I, Ch. 2, Stock No. 008-022-00113-1, 
U.S. Government Printing Office, Washington, D.C., 1977. 


Operating Limits and Warnings 
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@ 
User Instructions | 
Program Title: 102A Linear Wave Approximation to Breaking Wave Height 
and Breaking Wave Angle 


INSTRUCTIONS 


ENTEG BOTIOM sLOoPE MCTANGENT) 
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STEP 


INPUT 
DATA/UNITS 
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BEAD Kd IN DisPLA Pai | 
PRESS R/S AND READ Nn Leh Goma 
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PRess B/s AnD Read Ks ear aa 
PRESS R/s AND eead Hy!’ (Ft oc wm erate 
PRESS R/S and Read Hp (ft o2& wm ere ee 
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@ 
User Instructions — -’ 
Program Title: 102A Linear Wave Approximation to Breaking Wave Height 
and Breaking Wave Angle 


CET 

ENTER UW = I8ft PRESS D 

ENTER T= Ssec PRESSE 
\ 


OTHERWISE 2 
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@ @ 
Program Listing 
KEY ENTRY KEY CODE COMMENTS STEP KEY ENTRY KEY CODE COMMENTS 
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KEY ENTRY KEY CODE COMMENTS STEP KEY ENTRY KEY CODE COMMENTS 
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R/S - INSERT IF PRINTER 15 UNAVAILABLE 
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STEP KEY ENTRY &EY CODE COMMENTS STEP KEY ENTRY KEY CODE COMMENTS 
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Program Description 


Program Title 103R Shallow-Water Wave Forecasting Equations (RPN logic) - 
Name T.L. Walton, Jr. Date 10/80 


Address Coastal Engineering Research Center 


; Kingman Buildin ip Code 
ey Fore Belvoir 2 Soe Virginia a? 22060 


Program Description, Equations, Variables, etc. 

This algorithm computes the wave height, H, wave period, T, and minimum dura 
tion, t, from input values of the water depth, d, fetch length, F, and 
adjusted windspeed, Ua, using equations (1), (2), and (3) of CETIN-I-6. Equa 
tions (1) and (2) are for constant water depth and unlimited wind duration 
and have been revised from equations (3-25) and (3-26) of the Shore Protec- 
tion Manual. Wave height and period in this algorithm are significant wave 
height and period. Algorithm uses English or metric system of units. 


REFERENCES 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 
Proteetton Manual, 3d ed., Vols. I, II, and III, Stock No. 008-022-00113-1 
U.S. Government Printing Office, Washington, D.C., 1977, 1,262 pp. 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, "Method 
for Determining Adjusted Windspeed, Ua, for Wave Forecasting," CETN-I-5, 
Fort Belvoir, Va., 1981. 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, "Revised 
Method for Wave Forecasting in Shallow Water," CEIN-I-6, Fort Belvoir, 


Va., 1981. 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, "Revised 
Method for Wave Forecasting in Deep Water," CEIN-I-7, Fort Belvoir, Va., 
1981. 
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Program Description 


Program Title 103A Shallow-Water Wave Forecasting Equations (AOS logic) 

Name W.A. Birkemeier Date 5/81 
Address CERC Field Research Facility 

City Kitty Hawk, StaeNorth Carolina ZipCode 27949 


Program Description, Equations, Variables, etc. 

This algorithm computes the wave height, H, wave period, T, and minimum 
duration, t, from input values of the water depth, d, fetch length, F, and 
adjusted windspeed, Ua, using equations (1), (2), and (3) of CETN-I-6. 
Equations (1) and (2) are for constant water depth and unlimited wind dura- 
tion and have been revised from equations (3-25) and (3-26) of the Shore 
Protection Manual. Wave height and period in this algorithm are significant 
wave height and period. Algorithm uses English or metric system of units. 


REFERENCES 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 
Protectton Manual, 3d ed., Vols. I, II, and III, Stock No. 008-022-00113-1, 
U.S. Government Printing Office, Washington, D.C., 1977, 1,262 pp. 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, ‘Method 
for Determining Adjusted Windspeed, Ua, for Wave Forecasting," CETN-I-5, 
Fort Belvoir, Va., 1981. 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, "Revised 
Method for Wave Forecasting in Shallow Water,'' CETN-I-6, Fort Belvoir, 
Valeri 9 Sile. 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, "Revised 
Method for Wave Forecasting in Deep Water," CETN-I-7, Fort Belvoir, 
Va., 1981. 
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Program Description 


104R Depth-Limited Design Breaking Wave Height at Structure 
Program Title (RPN logic) 

Name Tate) Walton Jiy. Date 10/80 
Address Coastal Engineering Research Center 

Kingman Building i Zip Code 
Fort Belvoir, Virginia 


ci 22060 


Program Description, Equations, Variables, etc. 
This algorithm computes the depth-limited breaking wave height at a struc- 
ture for design purposes. It can be used in lieu of Figure 7-4 of the Shore 
Protection Manual. The equation for the curves in Figure 7-4 is not given 
in the SPM but can be found by simultaneous solution of SPM equations (2-91), 
(2-92), (2-93), (7-3), and (7-4). Input is wave period, T, and water depth 
at the structure toe, dg. The development of the equation is derived on the 
attached solution sheet. Algorithm uses English or metric system of units. 


REFERENCE 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 
Proteetton Manual, 3d ed., Vols. I and II, Chs, 2 and 7, Stock No. 008- 
022-00113-1, U.S. Government Printing Office, Washington, D.C., 1977. 


Operating Limits and Warnings 
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SOLUTION SHEET FOR PROGRAM 104R 


The following equations are given in the Shore Protection Manual: 
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Equation (7-4) can be rewritten in dimensionless form as: 


a 


A d 


Hp = = 


(b = al)! = mt. 
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Hp = Hp/gT2 and dg = dg/gT? 


The above equation can then be solved via the quadratic formula for Hy, in terms 


of ds, tp, m, a, and b where the positive root provides useful results. 


A ~~ A aA 1 =] 
Hp = { (mt pb - adg - 1) + [(mt pb-adg-1)? + 4amty bd | /2\. (2amtp) 


This is the equation used in the program for design breaking wave height. 
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Program Description 


Program Title ans Poison Design Breaking Wave Height at Structure 
AOS logic 
W.A. Birkemeier Date 5/81 


CERC Field Research Facility 
City Kitty Hawk, State North Carolina ZipCode 27949 


Name 
Address 


Program Description, Equations, Variables, etc. 


This algorithm computes the depth-limited breaking wave height at a struc- 
ture for design purposes. It can be used in lieu of Figure 7-4 of the Shore 
Protection Manual. The equation for the curves in Figure 7-4 is not given 
in the SPM but can be found by simultaneous solution of SPM equations 
(2-91), (2-92), (2-93), (7-3), and (7-4). Input is wave period, T, and 
water depth at the structure toe, dg. The development of the equation is 
derived on the solution sheet included with program 104R. Algorithm uses: 
English or metric system of units. 


REFERENCE 
U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 


Proteetton Manual, 3d ed., Vols. I and II, Chs. 2 and 7, Stock No. 008- 
022-00113-1, U.S. Government Printing Office, Washington, D.C., 1977. 
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Program Title: 104A Depth-Limited Design Breaking Wave Height 
at structure 


INPUT OUTPUT 
DELI SUS DATA/UNITS Oo Riakevss fy DATAJUNITS 
ENTER ScTICM SicrE m(tangent> m 
2 ENTER depth AT stRucTURE ds(fFt ce METERS) ds 


RE ee 
for (4) 
PRESS E' FOR METRIC UNITS for Cm 


NOTE4 WITH PRINTER M,ds, 7, Hb ARE 
PRINTED AND LABELED 


Exvampies ds=loh, T=!10sec, m= 0-10 


wn 
+ 
m 
a) 


L 
— 


ENTER m= G.-) PRESS A 
ENTER ds=10 f+ PREss B 10 Ft 
ENTER T=)0 SEC PRESSE IO SEC 


READ Ho = 17,98 IN. DISPLAY 
PRINTOUT FROM EXAMPLE : 


JUULUCUUUCUUUUUU See ee 
JUCUUUUBUUUCUUUUUU UU UU UCC 


lilies ee be Reed Jeu 


55 


STEP 


KEY ENTRY 


Program Listing 

a C2 a aaa Jew gay 
Cc « 

KEY CODE COMMENTS STEP KEY ENTRY KEY CODE COMMENTS 


v 


J 
r 


a tet Lene is WS 
1 nr Pe LEL 
$25 S00 MR, Le 

1 a yt = 
rere Es fo n> 3g Teens cove 
See = yA" PRINT CODE sal ! 
mith ry Sih OP og 
B39 OF 2 ales ines 
O¢ Od +* Sly (i 
45 RCL + OS) BE) es 
oo Oo ow Mes A OF m 
69 OF & Pent Bottom OSE ME eo rN RY v 
OSs 08 w« sure m G65 45 ROL 
31 Fes eS (Re ee 
Re (Leys Der) Bae 
ae cle 43° REL 
4 S70 ips Muss 
o nea d. ine 5 = 
o 1 Rah Wee: $3 ECL 
o 6 a DS PENT CODE a os 
0 hectw eo GS Se wee 
0 is or4 4 = 
5 FP 


* = 
Oo & Oi: 42 3570 due Ry 
o tg « nrée oO O41 “gt? 
ECL « Dee ce Rie 
a OF * OF 0 oo 
ee) (lP 7 i 5 EA 
Gen eet poe ds o i 
“ES O a 


LBL E' ALLows METRIC 


Wye) 


Ho Be ee ee ON ee PN 


OC A A a A ee OR 
mmo am 


I 
i E? DATA ENTRY a3 5 ten 
42 570 T>Ra Me4 22 THY 
os og eS) Semele 
ie) i S34 o+¢e- 
re S 9= 4.81 a Ps + 
ms -%5 0 o 1 
1 1 iM 5 = 
nies =) g>Re EM) ts 8 
ne of MSL 4 
me ¥ = mae 0 3 
woe. i * Sey) A 
E'S OP r Dee OR oF 
mT. mo (3c mL (> - 
Dees Aas ol® ee 43.45(I-e > 
63 OF *“ Pent “M* 10 WW ae Neat] = 
Os Wl Oe erie ere cen. WS eee ee eaAE PRz 
i3 Dg! Se (ie ns 
fia |Elete Was ae Ree 
I [S E - ALLOWS ENGLISH 100 00) 6D 
42 270 DATA ENTRY SSL es 
hey (He) TR Moe Me a 
(oie Ine eS) 
ire g=32.2 Ft/sec* EC SI a 
ENB} 5 peel mE O15 = 
De 2 Wiis: BIE 
eS > Rg te 


104A-3 


*-DELETE IF PRINTER IS NOT AWWAILABLE 


56 


STEP EY ENTRY KEY CODE COMMENTS STEP KEV ENTRY KEY CODE SUE COMMENTS 


= 
an 
if 
oe ¢ 


a 
eat ~)e 


ea 0) 
one 
(ee Cann iy 
el 
pau 


A 
4am % ds 


Ca) 


2] 
Niet eS 


7 [— 


(ee Oem Oe me 


, Y. 
(cR) +4 am tp) 


Goto fo 4e() c 


1.56 = 
3 a ve 74.5 ™) 


Ma OW Oe ed 


Da a se) 


5 tp 0) 6 


et eo 


» 

SRI GO a San 
i= 

ho 


APD ee Te oo eT A fe fl 
> wm 


om a ae) 


moe Be Be Pe OR ats ON ad ee OTR ee et 


fo 
1 
DOE ims 


symm ee ed ed ed Te TTT 


= 


LOOVEe WMO OOK OE 


Ca 8 6 


Cam 
ch 


1 CAN Lg ee a 


ag ve fe oe 


1m Oo 
» 
er 


0 {2 07 
OR Oe OOOO ee Oh oO no ho eee 
™ 


age my Bx Te Co me ae 
oO 


Pans 


‘\B peImy CODE 
m Tp = (4-4,.25m)M : 


oO 


>Re 


BO OO ho oo ho Po fe Pee he ee Pee Bo 


ho eo tt ee pe 


2 
3 
a) 
1 
2 
3 
4 
5 
6 
r 
& 
) 
a 


He 


PRINT Hy 


PSDP rr Oar a a Se SO 


Woh O Poe OMAN ePohe Com) Oe oof 


Oy OWS OY Oe Oe OO OO 


op rar) 
es | 


“, 


ae ee SL 


Aen 
a ¢ 


Dal a 
eee i 
(=) fu Ts 


a 
CT. i= 


104A-4 


Ye- DELETE |F PRINTER 15 NOT AVAILAGLE 


57 


Program Description 


105R Wave Transmission - Fuchs' Equation (RPN logic) 
T.L- Walton, Jr. Date 10/80 
Coastal Engineering Research Center 
Kingman Building 5 
Zip Code 
Fort Belvoir irgini fh 22060 


Program Title 
Name 
Address 

City 


Program Description, Equations, Variables, etc. 

This algorithm computes wavelength, L, in water depth, d, given the wave 
period, T. The program then computes wave transmission over a thin verti- 
cal barrier in water depth, d, using Fuchs" equation: 


where Ht is the transmitted wave height, Hy the incident wave height, and 
h the height of barrier. Note that this equation cannot be used when wave! 
transmission is by overtopping of a structure. Algorithm uses English or 


metric system of units. 


REFERENCE 


U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 
Protection Manual, 3d ed., Vol. II, Ch. 7, Stock No. 008-022-00113-1, 
U.S. Government Printing Office, Washington, D.C., 1977, p. /-62. 
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INSTRUCTIONS 
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Program Description 


Program Title = 105A Wave Transmission - Fuchs' Equation (AOS logic) 
Name W.A. Birkemeier Date 5/81 
Address CERC Field Research Facility 


City Kitty Hawk, State North Carolina ZipCode 27949 


Program Description, Equations, Variables, ete. 


This algorithm computes wavelength, L, in water depth, d, given the wave 
period, T. The program then computes wave transmission over a thin ver- 
tical barrier in water depth, d, using Fuchs' equation: 


where Ht is the transmitted wave height, Hy the incident wave height, and 
h the height of barrier. Note that this equation cannot be used when wave 
transmission is by overtopping of a structure. Algorithm uses English or 
metric system of units. 


REFERENCE 
U.S. ARMY, CORPS OF ENGINEERS, COASTAL ENGINEERING RESEARCH CENTER, Shore 


Protection Manual, 3d ed., Vol. II, Ch. 7, Stock No. 008-022-00113-1, 
U.S. Government Printing Office, Washington, D.C., 1977, p. 7-62. 


Operating Limits and Warnings 


105A-1 


6| 


User Instructions 


Program Title: 


INPUT 
INSTRUCTIONS | oaraunirs | 
ENTER WATER pe PTH (d) im FEET oR METERS | d( f+) ch 


; ENTER Sitt HEIGUT(h) IN FEET GR METERS h( Fpor 


T (sec 
AND LCABELED, Ly AND L ARE CCMPUTED IN 

READ Log 

PRESS R/S READL 


Seni Gee BeAS TORO i ia 
To REPEAT FoR NEWT GO TO STEP ener 


ve - 
0 


PRI 


d= 15 £ h_=\0F+ , T=IO sEC 


NTEX d=I5 pressA 
NTEP h=10 press B 


READ Ly = 512,48 
PRESS s REQuieeED READ L2 713.0 
PRES REQuIGED) PEAD U4 f 


JUUUUULUBLIULU Ue EOC OCUJOUO0000 


JUCUBUCU UC UCU UCU UUEEUE EO 


105A Wave Transmission - Fuchs' Equation 


OUTPUT 
OATA/UNITS 


d 


h 


Lo 
L 
tte Ji 


Pe) 
ees 28 
ee 8 7) 
ra Se 
pane 07 ©) 
2 oy 
ise 
gael 
ae S| 
a 


105A-2 


Program Listing 


STEP KEY ENTRY KEY CODE COMMENTS STEP KEV ENTRY KEY CODE COMMENTS 


Gi tisk i 
Lh AF wen | c 
See &) Ol awa Ro IU RSamean Re: 
oo oo 42 270 g7Re 
ee 3x Oe és 
ot 1 We Gia. (Usle 
Ge 6 ©  "D” pRit coe ae 
Ba OF 0 Bahn FF “1” PRINT CODE 
Ha et & meas 13 Sing 
Hoe) Pe 2 Med 6 a 
Hit Ca HES ol + 
mit 8S oq Hee 42 RL 
mis a5 = fer oO nF? 
His GHe [erste 217d 7Ro GES 6 
oa a0 ie a 
eo OF *  pewt d 
Lie E oy 
Si ees $2 RCL 
7e LBL Ge OG 
te. B da = 
42 S70 h7Rs ae 2 
os 08 ei 
ges VO He Sh 4 canny Gees 
Oe. = * “sy PRINT CODE Wy ° & 
a Me 2°92 (Stik HEIGHT) ey Sey 
a fl cc se He 2. & 
fj tae (DP te es GP ss 
r tre lid * aces dd * 
& ele Rel Be 
2 eS eS 5 = 
i Ge PS f 42 570 Le >Ry 
i G6 6 « PRINT h a) (a 
21 Res ee) ip 
ae (el steps 53 10 49 4 06 * ee re 
lm (e2 ALLCW Me TIC Fe LBL 
fer Scan INPUT ZOUTPUT hae 
ie ie T> Re ey a 
oa 4 Be ts 
ga g= 9.81 i 42 RCL 
m3o Oe 8 i mi oi 
Mah mii i Pe dell 
O41 42 570 g7?Re i; iz nt 
4g 0B UG <2 ST petow a? 
es- TS * p*® PRINT CODE ice ey ices 
es im a * cy (te R 
Des ES DF  S PRINT “M™ TO Sea. Be 
ne6 m4 fa * IWDICATE METRIC 43 REL 
were test (Be Ta 
nse  fS Se =o ao = 
mes «se [Re 42570 Znd >R4 
nso Te LEBEL SIE PS 50 1054 mies me howd 
oo iS E ALLOW ENGLISH ie I 
nes 32 Sag iNnPUT ZOUTPUT 
OSS ro MF TORY 
5 rie 3 


%- DELETE IF PRIMER 15 NOT AVAILABLE 10oe 
R/S- INSERT IF PRINTER 1S NOT AVAILAQLE 


63 
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Program Description 
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Address 
City Zip Code 


Program Description, Equations, Variables, etc. 
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